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ABSTRACT 

We present HST/NICMOS images at 0.2" resolution (Nyquist sampled in 
the central 19" ) of the HI Pa emission line in a 70" (a) x 90" (5) region of the 
Galactic center centered on the non-thermal radio source Sgr A*. The majority 
of the emission arises from ionized gas in the mini-spiral in the central parsec. 
Pa emission is also seen from 26 stellar sources, presumably early-type stars 
with mass-loss winds. The new data reveal significant small-scale structure 
(< 1" ~ 1.2 X 10^^ cm) in the ionized gas of the mini-spiral; low surface 
brightness emission features are also seen for the first time. The ratio of 
observed Pa emission to 6-cm radio continuum emission is used to derive an 
extinction map for the ionized gas over the central parsec. Ay varies from 20 
to 50 mag with a median value of Ay = 31.1 mag, in excellent agreement with 
earlier estimates derived for the stellar sources. Extinctions, independently 
derived using H92a recombination line data, were in excellent agreement with 
those derived from the 6-cm continuum. A broad minimum in the extinction 
extends ~ 30" NE-SW over the area of the IRS 16 cluster, including the area 
of peak flux in Pa and radio continuum. Large increases in extinction are 
seen along the periphery of the ionized gas, particularly in the direction of the 
western arm, suggesting that the ionized gas is partially extincted by dust in the 
molecular clouds at the outside of the ionized region; however the HII cannot be 
entirely behind the molecular clouds since the extinctions would then be much 
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greater. The minimum Lyman continuum emission rate within 40" radius of Sgr 
A* is 3.9x10^° sec"^ or L^^c - 2.7 x 10^ Lq, half of this within 20" radius. 
The small-scale, filamentary structures in the ionized gas have a free thermal 
expansion time of only ~ 3000 yrs; either magnetic fields or mass-loss winds 
from the hot emission line stars may contain the ionized filaments. (Eight new 
emission-line stars were also detected in Pa.) 

For both the ionized gas and stellar continuum, the centroids of the emission 
(corrected for extinction) remain within ~ ± 1" from a radius of 2" out to 
40"This provides further evidence that Sgr A* (the putative central black hole) 
is indeed at or extremely close (< 0.04 pc or 1.2 x 10^^ cm) to the center of the 
Galactic nucleus stellar distribution and presumably the dynamical center of 
the Galaxy. The 1.9;um surface brightness increases inwards to 0.9" radius and 
exhibits a decrease or leveling off closer to Sgr A*, possibly indicating the core 
radius of the central stellar distribution or depletion of the late-type stars by 
stellar collisions near the central black hole. 

Subject headings: Galaxy: center — dust, extinction — infrared: ISM — HII 
regions 

1. Introduction 

Within the central few parsecs of our Galaxy, radio, infrared and x-ray observations 
have revealed a massive black hole lying within a region of ionized 'spiral' arms, encompassed 
by a ring of dense molecular clouds. The black hole is identified with the non-thermal radio 
source Sgr A*, and it is estimated to be 2.6x10^ M© from the motions of nearby stars and 
the neighboring ionized gas (Eckart & Genzcl 1996, Gcnzcl et al. 2000, Ghez et al. 2000). 
At a projected radius of approximately 2 pc from Sgr A*, lies the circumnuclear disk/ring 
(GND) of dense molecular gas and dust clouds (Genzel et al. 1985, Giisten et al. 1987), 
and within this ring, ionized gas streams extend down to a few arcsec from Sgr A* (Lo & 
Claussen 1983; Lacy et al. 1991). The molecular cloud ring can be modelled as a torus (0.5 
pc thick) inclined at 70° to the line of sight (Jackson et al. 1993). Although clearly much 
less energetic and massive than the luminous AGN seen in many galactic nuclei, this region 
probably provides our best opportunity to observe in detail the processes associated with 
buildup of massive black holes and their accretion processes. 

Previous studies of the ionized gas in the inner few parsecs of the Galaxy, with an 
emphasis on the radial velocity kinematics have included observations of [Nell] (Lacy, 
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Achtermann, & Serabyn 1991), Br7 ( Burton & Allen 1992 and Herbst, Beckwith, & Shure 

1993), and H92q; (Roberts & Goss 1993). Roberts & Goss published a ratio map of their 
H92 a radio observations to Br7 from Burton & Allen (1992) which gave a qualitative 
representation of the extinction in the inner few parsecs, however, they did not calculate 
the extinction in Ky. 

We have recently undertaken two investigations aimed at providing the highest 
resolution imaging of both the ionized and molecular gas in the circumnuclear disk. Most 
of the existing observations of the ionized and dense molecular gas are at > 1 and > 
8 " resolution, respectively, and thus do not clearly resolve the internal structure of the 
ionized gas. The molecular gas clouds were imaged in HCN at 3" resolution by Christopher 
et al. (2003) using the Owens Valley mm-Array; here we present HST-NICMOS imaging 
at 0.2" resolution of the 1.87//m Pa (HI 4-3) recombination emission line. The Pa line 
has 12 times greater intrinsic flux than the 2.16//m Br7 line, and thus it provides excellent 
sensitivity when observed with the NICMOS camera on HST (despite the extra magnitude 
of extinction at the Galactic center at 1.87/xm compared with 2.16/^m). Comparison of the 
Pa emission with radio continuum emission from the same ionized gas is used to derive a 
map of the foreground extinction; previous estimates of the extinction have been mostly 
derived from the individual stellar colors; the derived extinctions are therefore not spatially 
continuous and their location along the line of sight is often uncertain. For both the HCN 
and Pa observations, multiple fields were observed in order to image the entire Galactic 
center region out to and including the CND. 

In §2, we describe the observations and reduction procedures. Three methods for 
removing the stellar continuum from the emission line images (scaled subtraction of the 
190N image from the 187N image, PSF fitting and subtraction for the bright stars and 
variable scale factor image subtraction) are discussed in Appendix 1. In §3, we derive 
the extinction distribution from the observed ratio of Pa to radio continuum and radio 
recombination lines and present an extinction-corrected Pa image. Estimates of the sizes 
and densities of the Hll structures are derived in §4, and the stellar sources with emission 
line excesses are discussed in §5. We then examine the observed and extinction-corrected 
diffuse light distribution of stars in §6 and the stellar density profile relative to Sgr A*. 
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2. Observations and Data Analysis 

2.1. NICMOS Pa Observations 

Images of the CND were obtained with the NICMOS camera on HST in the 
narrow-band 187N (Pa) and 190N filters as part of the NICMOS GTO program. A mosaic 
image, extending out to 50" radius, was constructed from four Camera 3 (NIC3) dithered 
fields and one Camera 2 (NIC2) dithered field (for higher resolution near the center). N1C3 
uses a 256x256 HgCdTe array with plate scale of 0.203" per pixel providing a 52.0" field of 
view; Camera 2 has 256x256 0.076" pixels and a 19.3" FOV (Thompson et al. 1998). The 
F187N and F190N filters with effective wavelengths of 1.87 and 1.90 /im were used to obtain 
on-line and off-line images. The F187N filter has a bandwidth of ~ 1% corresponding to 
over 3000km s~^; it therefore includes all velocities of the Galactic center ionized gas (Lacy 
et al. 1991). The filter also includes a HeI(4-3) doublet fine (at 1.8691 and 1.8702 //m) ; 
for the extended emission. Pa is expected to completely dominate (by more than a factor 
of 30), but on some emission lines stars (Paumard et al. 2001) the Hel fine could be a 
significant contribution. 

For each of the field centers (listed in Table 1), 4 dithers were taken spaced by l.l"with 
non-destructive reads at each position (MULTIACCUM). The total integration times for 
each filter are listed in Table 1. The diffraction limited resolution of HST is 0.19" at 1.87yum; 
the N1C2 images are therefore sampled at better than the Nyquist frequency but the NIC3 
images are undersampled by the Nyquist criterium. Although not completely recovering 
the Nyquist sampling, the half-pixel dithers do help to reduce the undersampling in NIC3. 

2.2. Data Reduction 

The data were calibrated using the CALNICA version 3.3 task (Bushousc & Stobic 
1998) in IRAF^/STSDAS with most reference files (static data quality, detector read noise, 
detector non-linearities files) from the Space Telescope Science Institute (STScl) NICMOS 
pipeline. The flat-field and dark frames were provided by the NICMOS Instrument 
Development Team (IDT). The dithered images were then shifted and coadded; for the 
NIC2 fields, using the the IDP3 routine from the NICMOS IDT software. We estimate 
the registration accuracy between the NIC2 fields to be better than 0.05 NIC2 pixel. The 
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registration accuracy for the NIC3 fields was less accurate due to the undersampling of the 
point spread function (PSF) and was of order 0.25 NIC3 pixel (0.05" ). The NIC3 data 
were dithered by half-pixel steps to oversample the original pixels by a factor of 2 in the 
final mosaic. The NIC2 and NIC3 mosaics were rotated to north up and east to the left 
using the data-header orientation angle and uniformly sampled to 0.0753" square pixels, 
corresponding to 0.0029 pc or 9.00 x 10^^ cm at a distance of 8.0 kpc. (The pixels are 
factors of 1.00906 (NIC2) and 1.005 (NIC3) larger in x than in y - this was corrected to 
make square pixels by interpolation). Where the NIC2 and NIC3 fields overlapped, offsets 
were determined from centroiding at least 20 common stars in the overlap regions. In the 
central area covered by NIC2 with better sampling, the final mosaic was made from just 
the NIC2 pixels. 

For the NIC3 images, we employ scale factors of 5.050x10"^ and 5.033x10"^ 
Jy (ADU/sec)-^ at 1.87 (Pa) and 1.90 /xm, respectively (Rieke et al. 2003). For the NIC2 
images, the scale factors are 4.086x10^^ and 4.330x10"^. The rms noise measured in the 
final NIC3 and NIC2 images in areas of very high extinction with relatively few stars is 
typically 3.9 and 2.15 /iJy per pixel, respectively. 

The mosaiced F187N and F190N images are shown in Fig. 1 with a logarithmic stretch 
on the grey-scale. In the F187N image diffuse emission (Pa) at the peak level of ~ 200 
/iJy per pixel may be seen; however, over 850 stars in the image are at a level above this 
brightness. Also apparent in these images arc several areas with far lower star counts 
than their surroundings - presumably due to small but very opaque dust clouds in the 
line of sight. Contours of HCN (1-0) emission at 3" resolution (Christopher et al. 2003) 
are superposed on the 1.87 /im image; several areas of decreased stellar density are well 
correlated with molecular clouds seen in HCN (see Fig. 1), implying that the dense dust 
clouds are located near the Galactic Center as pointed out in an early analysis of this work 
in Stolovy et al. 2001. 

2.3. Continuum Subtraction and Emission Line Flux Calibration 

In order to image the Pa line emission, it is necessary to remove the much stronger 
stellar point sources. The removal of the stellar continuum from the F187N images is 
difficult since the stellar density is high, the intrinsic colors of the stars vary and most 
importantly, the reddening is highly variable across the field. We experimented with three 
techniques for continuum subtraction (subtraction of fixed-scaled version of the F190N 
image, point source removal of the stars, and subtraction of a variable-scaled version of 
F190N) and the results are discussed in the Appendix to this article. The final results 
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presented below are based on the last technique. A list of emission-line stars is provided in 
one of the sections below. 

To convert the flux densities in the F187N continuum-subtracted image into an 
emission line flux, we adopted an effective FWHM for the filter of 190A or /S.v = 1.63 x 
10^2 Hz. 



3. Pa Images 

Fig. 2 shows the Pa emission for the entire region imaged. The reference position of 
Sgr A* ( asooo = nHh"^ 40^0383 ± 0.0007, ^2000 = -29°00'28.069 ± 0.014"; Yusef-Zadeh 
Chaote & Cotton 1999) was located on the NICMOS image using the offsets from maser 
source IRS 7 given in Menten et al. (1997). The relative offsets of all other maser sources 
specified by Menten et al. (1997) (IRS IINE, IRS 7, IRS 15 NE, IRS 10 EE, IRS 16 NW, 
IRS 16 NE and IRS 15NE) agreed to within one 0.0753"pixel in the final mosaic image. 
Subtraction of a scaled F190N mosaic from the F187N mosaic should in principle yield 
an image of the Paa emission, but, as shown in Fig. 13, strong positive and negative 
stellar residuals (due to varying intrinsic colors, varying extinction, and photospheric Pa 
emission and absorption) are apparent. We instead removed the continuum with a scale 
factor that varied spatially across the field to compensate for these issues as described in 
the Appendix. This variable spatial scale factor continuum removal method was employed 
in the subsequent Pa figures and anlaysis presented here. 

The most prominent large-scale emission features in Fig. 2 correspond to : the 
East- West 'bar' structure approximately 7" south of Sgr A*, the Northern arm which 
extends from -2" to over 25" north of Sgr A*, and the western arm approximately 14" west 
of Sgr A*. With the high resolution and sensitivity of the Pa images, these large-scale 
features break up into numerous filaments/arcs which are presumably density enhancements 
in the ionized gas or individual ionization fronts at the edges of neutral clouds. Smaller 
scale emission peaks arc seen on the 26 emission line stars. The 'mini-cavity' appears as a 
clearly evacuated spherical hole in the emission of diameter 2.7" at 1.7" W, 2.8" S of Sgr 
A* (Yusef-Zadeh, Morris, & Ekers 1990) and is as shown in Fig. 2. A weak Pa emission line 
star appears at the center of this cavity (Stolovy et al. 1999). Another weak emission line 
star is seen at the center of the "eastern cavity" (4.6" E, 3.4" S ). 
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3.1. Extinctions from Pa, Radio Continuum and Radio Recombination Lines 

The free-free and recombination line emission at radio wavelengths suffer negligible 
extinction. The observed ratios of Pa-to-radio free-free (and recombination line flux) can 
therefore be used to estimate the dust extinction - both the line and free-free emission 
depend of the line-of-sight integrated emission measure, but only the Pa line emission is 
significantly extincted by dust. 

For case B recombination, the intrinsic Pa flux in a pixel is given by 



where Te, Ug and Up are the temperature and volume densities of electron and protons, 
1 is the path length in the ionized gas, apj,j is the projected area of a pixel and d is the 
distance (Table 4.4 in Osterbrock 1989; the Pa emissivity used to derive the constant in 
Eq. 1 is appropriate to ng = 10^ cm~^ but the dependence on Ue is very weak.) In Eq. 1 
and the following, we will adopt Tg = 6000 K (see Morris & Yusef-Zadeh 1989; Simpson 
et al. 1997). Similarly, the radio free-free continuum flux density from a pixel is given by 

where v is the frequency of the radio emission. Thus the intrinsic Pa/free-free ratio is : 



( EfA ^ 1.54 X 10-^=^ ( _I^y'''' ( ' erg s-^ cm-^ mJy-^ (3) 

Similar to the radio free-free continuum emission, radio recombination line emission 
also traces the same ionized gas as Pa but suffers negligible extinction, providing an 
alternative means of estimating the extinction. Roberts & Goss (1993) observed the H92a 
recombination hne from Sgr A West. If the line is emitted under LTE conditions (Roberts 
& Goss 1993), then the line flux integrated over the velocity width of emission per pixel is 
given by 



5...A. = 1.29xlO-(^)- 3,3;,/':;';'° mJykn.s- (4) 
using Eq. 13.28 from Rohlfs & Wilson (1996). The intrinsic Pa/H92a ratio is therefore : 



Sm2a^v J intrinsic VeOOOXy \%.?,Ghz 



(5) 



The observed flux ratio can be less than that in Eq. 3 and 5 due to extinction at A = 
1.87^m (Pa). Rieke (1999) derived Ai.go^m/Ay = 0.176 and A2.22,im/Ay = 0.105 from 
NICMOS observations of stars in the Galactic center; interpolating between these values 
with a A~^'^ power-law yields Ay = 7.24 x Kpa- (This is nearly equal to the value obtained 
from Rieke & Lebofsky 1985.) We therefore use the relation 



to estimate the visual extinction where SffH92a is the free- free or H92q; flux. 

To estimate the extinction of the Pa emission, we use the A = 6 cm radio map of 
Yusef-Zadeh & Wardle (1993). This VLA radio image at 0.40" x 0.67" resolution has 

good sampling at low spatial frequencies; the image therefore contains nearly all the 
emission on scales corresponding to the CND Pa emission. The Pa image was smoothed 
to the resolution of the radio image before computing the observed ratios of Pa to radio 
continuum. 

Fig. 3 shows the extinction distribution obtained from the observed ratios of Pa to 
radio continuum using Eq. 6 and a histogram of the derived extinctions is shown in Fig. 4. 
The extinctions range from Ay ~ 20 to 50 mag and the median pixel extinction is Ay = 31.1 
mag. The median extinction is in good agreement with the median derived from the colors 
of the Galactic center stars (Blum et al. 1996, Cotcra et al. 2000). However, the extinction 
exhibits large-scale gradients across the region (see Fig. 3) and low extinctions (20-25 mag) 
are seen in front of the East- West bar and the Northern arm and in the direction of the IRS 
16 cluster. The generally low extinctions seen along the NE-SW swath south of Sgr A* is 
consistent with the low extinction values derived for the stellar continuum in IRS 16 (Blum 
et al. 1996). High extinctions (35-45 mag) arc seen in the Western arm and in general along 
the outskirts of the Pa emission. This may be due to the increased extinction associated 
with dust in the molecular gas along the periphery of the ionized region. The apparently 
low extinctions toward the Pa emission line stars in the central parsec are contaminated by 
photospheric excess Pa and should be disregarded in interpreting the line-of-sight extinction 
in the diffuse gas. The apparently high extinctions shown on Sgr A* and a few arcsec to 
the north are spurious, due to the presense of non-thermal radio emission there; the Ay 
derived from the H92a hue shows no anomaly on Sgr A* (see below) . 
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The extinction derived from the 6-cm radio continuum image suffers from two major 
uncertainties : non-thermal emission and the removal of spatially smooth, low-level emission 
in the VLA maps. The non-thermal contribution is significant on Sgr A* and a small 
contribution could arise from Sgr A East (Pedlar et al. 1989). For this reason we also 
derived the extinction using the H92q; radio recombination line image obtained by Roberts 
& Goss (1993). This H92a data, smoothed to 2" resolution, has a considerably lower 
SNR (and spatial resolution) than the 6-cm image, but it is completely immune from 
problems of non-thermal emission contamination by Sgr A*. Fig. 5 shows the extinction 
distribution obtained from the observed ratios of Pa to H92q; emission using Eq. 6. The 
areas with sufficient signal-to-noise ratio exhibit extinction values in close agreement with 
those calculated from the 6-cm data (see Fig. 3). Using the H92q; data, the median pixel 
extinction is Ky = 29.2 mag, compared with 31.1 mag from the 6-cm analysis (see above). 
In Fig. 6 we show a histogram comparing the extinction extimates from H92q; and 6-cm 
radio continuum on a pixel by pixel basis, illustrating the general consistency of the two 
approaches. Because of the high sensitivity, resolution and spatial coverage, we will use the 
extinction derived from 6-cm radio continuum in the discussion below, but the agreement 
between the two methods reassures us that there are not large systematic errors. In 
addition, the fact that both determinations exhibit a minimum in the ionized bar region 
and Northern arm argues strongly that this minimum in the extinction is real. Lastly, we 
note that we also used a 2-cm radio map (Yusef-Zadeh, Morris, & Ekers 1990) for a parallel 
extinction analysis and derived very similar estimates but this data has lower SNR than the 
6-cm data. 



3.2. Extinction-corrected Pa and Lyman Continuum Luminosity 

The distribution of Pa emission corrected for extinction at 1.87//m is shown in 
Fig. 7. Although the original radio continuum image is also, of course, extinction- free, 
the extinction-corrected Pa image shown in Fig. 7 has the advantage of higher spatial 
resolution and signal-to-noise ratio. Under the assumption that the extinction does not 
vary on scales smaller than the radio image resolution (0.40" x 0.67"), this image should 
provide the most accurate representation of the diffuse HII distribution. A histogram of the 
derived extinction-corrected Pa fluxes (per pixel) is shown in Fig. 8. 

The total fiux (within 40" radius of Sgr A*) in the extinction-corrected Pa image is 
2.4x10"^ ergs cm^^ sec~^; the luminosity is therefore 1.8x10"^^ ergs sec^^ or 1.7x10'^^ Pa 
photons per sec. For Case B recombination at 6000 K (Roberts & Goss 1993; Roberts, 
Goss, & Yusef-Zadeh 1996), aB{H)/aB{Pa) = 22.7 (Osterbrock 1989) and the total Lyman 
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continuum (LyC) emission rate required for the Pa emission is 3.94x10^'^ sec~^ or 2.70x10^ 
Lq (assuming 1.2 Rydbergs per Lyman continuum photon). Within 20" radius of Sgr A*, 
the Pa flux is 50% of that given above; thus 50% of the ionizing photons are absorbed 
within 20" or 0.8 pc. All of the emission line stars are within this region. The derived LyC 
production rates arc minimum estimates since we assume Case B recombination with no 
Lyman continuum escaping the region. Given the apparent geometry of a tilted ring in 
the neutral gas, it is, in fact, very likely that the actual production rate could be several 
times higher. The 8 emission line stars studied by Najarro et al. (1997) have a combined 
Lyman continuum emission rate of 3x10^*^ sec~^ and thus they can probably account for the 
ionization without any extra input from Sgr A*. Sgr A* has not yet been unambiguously 
identified as an infrared point source in either the continuum or lines. Extended emission 
is seen in our images at the location of Sgr A* but no discrete source is seen. The Lyman 
continuum emission rate obtained from Pa is consistent with that derived directly from the 
radio continuum (~ 3 x 10^° sec~^, Genzel, HoUenbach, & Townes 1994). 

4. Physical Conditions in the Ionized Gas 

4.1. Sizes of Ionized Structures 

The structures seen in the ionized gas are well resolved in the NICMOS Pa images. 
The emission filaments in the Northern arm and bar have typical FWHM ~ 0.3 - 1" (or 
0.4 - 1.2x10^^ cm~^) in their small dimension and up to 10 arcsec in their long dimension. 
The highly elongated, fine structure in the ionized gas is strongly suggestive that the 
structures are often ionization fronts on the edge of more extensive neutral gas clouds or 
arms. The ionized shell at the 'mini-cavity' has a diameter of 2.7" or 3.2 xlO"*^^ cm. Stolovy 
et al. (1999) discuss the structure of the 'mini-cavity' in Pa and [Fell] emission lines and 
identify the central star which they propose to be the source of the cavity. The age of the 
minicavity is only 500 yrs for assumed constant expansion velocities of 100 km s^^. It is 
possibly even a younger structure; Zhao & Goss (1998) estimated a current expansion rate 
of 200 km s~^ from radio proper motion measurements. 

4.2. Electron Densities (ug) 

Since the extinction-corrected Pa emission is proportional to the volume-integrated 
emission measure, the electron density in the ionized gas can be obtained from the 
extinction-corrected Pa image. Inverting Eq. 1, we find : 
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<ne> = 1.21 X 10^ F'Jl,, r°ooo T, ^ cm-' (7) 

where Fp^-u is the hne flux in 10~^^erg cm~^, Tgooo is the temperature normahzed to 
6000 K (Roberts & Goss 1993, Roberts, Goss, & Yusef-Zadeh 1996), Dg in the distance 
to the Galactic center (8.0 kpc), I17 is the linc-of-sight pathlcngth (normalized to 10^''' 
cm corresponding to 0.83" ), and a^j^ is the projected area of the pixel normalized to 
0.0753x0.0753" . The peak extinction corrected Pa fluxes along the Northern arm ridge 
are 4.3x10^^^ erg s~^ cm~^ per pixel. Adopting In = 1, the electron density in these peaks 
is Ue = 2.5x10"^ cm~^ and the mean along the Northern arm is n^, ~ 10'^ cm~^. These 
estimates are in good agreement with those originaUy derived by Lo & Claussen (1983) 
from the radio free-free emission. 



4.3. Relative Geometry of Ionized and Moleculair Features 

Some of the large-scale ionized structures can be ascribed to photoionization of the 
exterior molecular gas by Lyman continuum from hot stars in the central star cluster. The 
clearest example of this is the Western arm. It extends over 70" (N-S) in both the ionized 
and molecular gas tracers and the Pa emission centroid is everywhere offset 5-10" inside 
the HCN (Fig. 7). In the western arm, the ionized and molecular gas velocities are also 
similar (Christopher et al. 2003). On the other hand, the ionized structures at smaller radii 
(e.g. the bar and lower Northern arm) do not have any obvious molecular counterparts. 
These may be structures which are almost entirely ionized, elongated dynamically and 
compressed laterally by winds (see below). The north end of the Northern arm in Pa 
apparently connects smoothly in space and velocity (Lacy, Achtermann, & Serabyn 1991) 
to a molecular feature extending south (Aa = 10" , A5 = 30" ). In this instance, possibly 
the neutral gas becomes ionized at ~ 18" as it approaches the central cluster. Jackson 
et al. (1993) have suggested that the Northern arm is the ionization front on the surface of 
an atomic cloud seen in the 63//m [01] line, which fills the interior of the CND but displays 
a twin peak morphology similar to that seen in the integrated 2.12/im H2 emission line 
(Yusef-Zadeh et al. 2001) associated with the molecular torus of the CND. Further evidence 
for a neutral component of the gas interior to the CND was shown by Stolovy (1997) from 
[Sill] 34.8 /im observations at 15 " resolution, which peaked near the "minicavity" possibly 
as a result of enhanced Si abundance due to grain destruction in shocks. 

Correlation (or lack of correlation) between the molecular gas mapped by the HCN 
contours in Fig. 3 and the derived extinctions in front of the ionized gas might be used 
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to infer the relative placement of the H2 and HII gas along the hne of sight. The strong 
increase in extinction at the edge of the Western arm suggests that some of the ionized gas 
is behind, or mixed with, the neutral gas and dust in the molecular Western arm. Increased 
extinction in front of the ionized gas is also apparent where two filaments of molecular gas 
(traced by the HCN contours in Fig. 3) extend southward into the northern boundary of 
the ionized gas (Aa = -5" , AS = 15" , see Fig. 3). 

However, Christopher et al. (2003) estimate volume densities n^^^ > 10^ cm^'^ for the 
HCN molecular clouds. Their typical sizes (~ 5" or 7x10^^ cm) then imply column densities 
~ 7 X 10^^ H2 cm~^ or Ay ~ 700 mag for a standard Galactic dust-to-gas ratio. Clearly, the 
extinctions in front of the the Western arm and northern Pa emission are not nearly this 
high; if they were, no emission would have been detected in the near-infrared. We therefore 
conclude that Pa emission on the west and northern border of the region (Act = -10" , A5 
= 10" ) must be on the front face of the molecular gas (unless the dust is severely depleted 
or the molecular gas is extremely clumpy on scales << 1" ). 



4.4. Tidal Disruption and Stability 

The extraordinary small scale structures now seen in both the molecular and ionized 
gas in the central parsecs, must be extremely short-lived and evolving rapidly. The 
dynamical time (R/Vrot) is only 1.3x10^ yrs at R = 1.5 pc with Yrat — HO km s"^ (Jackson 
et al. 1993). If some of the ionized gas is infalling at ~l/2 of the rotational velocity as 
suggested by their shapes (for example the Northern arm as modelled by VoUmer & Duschl 
2001 and Sanders 1998), then the infall times are only a few 10"^ yrs. 

Similar short timescales are derived from the small scale structure in the Pa 
emission. Along the Northern arm, the width of the bright emission ridge is just 10^^ cm 
(corresponding to FWHM = 0.83" ); if this gas freely expands at the 10 km s^^ sound 
speed of the HII, the time required to double its width is just 3x10'^ yrs. Thus if the 
ionized gas is to maintain its high surface brightness over a dynamical or infall time, the 
observed emission ridges must be ionization fronts on the edge of dense, neutral material 
or the ionized gas must be compressed by a much lower density and higher temperature 
confining medium (e.g. a wind from Sgr A* or the mass-loss stars - see below) or magnetic 
fields. Most likely, both scenarios are occurring. The ionized gas on the Western arm is 
well-correlated both spatially and kinematically with the dense molecular gas seen in HCN 
emission (Christopher et al. 2003) and thus is likely to be the ionized inner face of the 
molecular arm. 
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A number of investigations starting with Quinn & Sussman (1985) have emphasized 
the importance of tidal disruption to the central parscc gas clouds. The Roche limit 
associated with the central 2.6x10^ Mq (i.e. the Sgr A* black hole; Genzel et al. 2000; 
Ghez et al. 2000) is n^ > 2.4 x IOVR^^ cm-^ (e.g. Sanders 1998). This is a factor of 10^ 
(!) greater than the density estimated above for the ionized gas in the Northern arm and 
is comparable to or greater than the densities derived by Christopher et al. (2003) for the 
HON clumps at ~ 2 pc radius. Inclusion of the mass associated with the stellar distribution 
doubles the Roche limit at the radius of the CND. 

The data clearly show the ionized gas in the Northern arm directly connecting to a N-S 
feature in the molecular gas. The kinematics of the [Nell] (Lacy, Achtermann, & Serabyn 
1991) and HCN emission are also consistent with this connection. Thus, the Northern arm 
is very likely an infalling streamer originating from the northern molecular gas (Jackson 
et al. 1993, Wright et al. 2000). The molecular gas might have been tidally sheared and 
ionized as it approaches Sgr A* and the hot central stars. 

The sharpness of the Northern arm in cross-section is illustrated in Fig. 9 where a 
1-pixel wide slice in the E-W direction is shown for the observed and extinction-corrected 
Pet. The width is less than 1" . Confinement of the ionized gas in the east-west direction 
requires a confining medium, magnetic field, or wind ram pressure. Paumard et al. (2001) 
show that there are two classes of emission line stars : those with broad emission lines 
{AV ~ 1000 km s^^) and those like the IRS 16 sources with widths of ~ 200 km s~^. The 
mass-loss stars in the central He cluster have a total M~ 3 x 10^^ Mq yr~^ and typical 
velocities v^ ~ 750 km s~^(Hall, Kleinmann, & Scoville 1982; Krabbe et al. 1991; Najarro 
et al. 1997). At a distance of 10" (0.4 pc) from the mass-loss stars, the wind will result 
in a maximum ram pressure pv'^r^ 9 x 10~^ dynes cm~^ (assuming no deceleration of the 
wind) . The thermal pressure (nkT) of the ionized gas in the Northern arm is approximately 
two orders of magnitude less (8xl0~^ dynes cm~^ assuming n = 10^ cm~^ and 6000 K). 
Thus, the mass-loss winds might reasonably contribute to the confinement of the ionized 
gas in the arm. It is also noteworthy that Baganoff et al. (2001) have recently detected 
X-ray emission from hot, ionized gas in the vicinity of Sgr A*. This gas could drive an 
outflowing wind, which might also compress the Hll gas in the Northern arm. Aitken et al. 
(1991) concluded from mid-infrared polarimetry that the magnetic field is aligned parallel 
to the flow of ionized material along the fllaments and that the magnetic fleld is therefore 
magnifled due to these shearing motions. They also concluded that the IRS 16 sources do 
not disturb the magnetic fleld and therefore they are hkely out of the plane of the Northern 
arm. 

Thus, it seems hkely that the Northern arm gas must be conflned in the El-W direction 
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by a wind or a confining medium from the general direction of Sgr A* (but not necessarily 
from the IRS 16 sources-the IRS13 wind sources and Sgr A*itself may contribute to the 
wind). The shape of the Northern arm, concave around Sgr A*, is also suggestive of pressure 
confinement on its west side. In the east direction (away from Sgr A*) the Pa emission in 
the Northern arm falls off more gradually and hence does not require confinement. 

5. Candidate Pa Emission Line Steirs 

We used the IDL 'starfindcr' (Diolaiti et al. 2000) program to derive the fluxes in each 
of the F187N and F190N filters by PSF matching for the NIC2 data only which led to 
satisfactory results in deriving stellar fluxes as compared to aperture photometry and scaled 
PSF subtraction. (The NIC3 portion of the mosaics, despite dithering and resampling of the 
pixels, did not have not sufficiently well defined PSF's to use PSF- matching techniques). We 
then correlated well-matched sources found in both filters and computed the F187N/F190N 
ratios. Although the large width (1%) of the NICMOS F187N filter makes it difficult to 
detect stars with relatively weak Pa line emission on top of the strong stellar continuum, 
we do list 26 stars 'candidate' emission line stars in Table 2 which were chosen based on 
derived F187N/F190N flux ratios > 1.1. The typical F187N/F190N ratio for stellar sources 
was found to be 0.91; a ratio of 1.1 is therefore a conservative indicator of significant 
(~ 20%) excess flux at 1.87;um. The high reddening of the Galactic center stars implies 
F187N/F190N < 1 for all stars. The fluxes in Table 2 have not been corrected for extinction. 

The last column in Table 2 lists the likely identiflcation of the stars with previously 
observed Hcl and HI emission line stars (Genzel et al. 1996; Najarro et al. 1997; Paumard 
et al. 2001). Several of these candidate Pa emission line stars are also strong Hel (2.058/im) 
emission line stars and the F187N flux may indeed be dominated by the Hel line within 
the bandpass rather than Pa. Some stars in Table 2 appear to be newly identified. (We 
note that some of the stellar positions listed here disagree with those given by Paumard 
et al. (2001) by as much as 0.8" ; we have discussed this discrepancy with Paumard 
et al. and their revised positions are in better agreement.) We are confident of our positions 
since the reference stars given by Mcntcn et al. (1997) agree to within 0.07" in the NIC2 
image. Several stars with high F187N/F190N ratios were removed from the candidate list 
after visual inspection determined that the excess was dominated by diffuse interstellar Pa 
emission along the line of sight. Weaker Pa emitters, such as IRS16CC, the 'minicavity 
star', and the 'eastern cavity star' do not appear in Table 2, as they do not obey the 
criterion set by F187N/F190N fiux ratios > 1.1. 'Starfinder' did not resolve the 2 emission 
line stars located near IRS13E, now called IRS13E2 and IRS13E4 (Maillard et al. 2003). 
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The fluxes for these stars were derived by a scaled subtraction of the PSF from each star in 
each filter. IRS15NE appears as 2 stars: a Pa absorption line star to the south coincident 
with the maser position of IRS 15NE and a star 0.5"to the north, which is an emission line 
star only imaged in NIC3. Aperture photometry was used to derive the Pa flux for the 
northern source, which we call IRS 15NNE. 

6. Centroid and Radial Distribution of Stellar Continuum 

It is of interest to determine if Sgr A* is located near the center of the nuclear 
star cluster and the ionized gas. Previous ground-based observations of the stellar light 
distributions have included the contribution from all stars, including the bright stars. Here 
we investigate the diffuse light surface brightness distribution by measuring the median 
pixel brightness as a function of radius from Sgr A* and the mean of pixels excluding the 
20% of the pixels with highest brightness. The very clean and stable PSF of NICMOS on 
HST is particularly advantageous for this analysis. We have also measured the centroids of 
the total Pa and 1.9//m stellar continuum flux. Although the starlight may not have the 
same effective foreground extinction distribution as the ionized gas, it is likely the stellar 
light corrected for the extinction derived for the ionized gas is a much better approximation 
to the intrinsic 1.9;um continuum than the observed 1.9/im continuum. For the extinction 
distribution we use that shown in Fig. 3 with the exception that within ± 0.5" of Sgr A*, 
the extinction was set to the average between 0.5 and 0.75"fadius which is within 10% 
of that determined from H92q;. The non-thermal flux of Sgr A* yields a spuriously high 
apparent extinction on Sgr A*. In computing the brightness distributions, we also exclude 
the very bright source IRS 7. At each radius in the distributions, we only include those 
pixels within the mosaic area with derived extinctions (and in the case of Pa only those 
pixels with detected Pa). 

The emission centroids for the total Pa and the 1.9//m flux were measured for apertures 
of radius 2, 5, 10, 15, 25, 30, 35, 40" from Sgr A*. These centroids are plotted as boxes in 
Fig. 10 on top of the 1.9//m image of the region within 2.5" of Sgr A*. For both the ionized 

gas and stellar continuum, the centroids remain within ± 2" over the entire range out to 
40" radius. This provides very strong evidence that Sgr A* (the putative central black hole) 
is indeed at or extremely close (< 0.04 pc or 1.2 x 10^^ cm) to the center of the Galactic 
nucleus stellar distribution and presumably the dynamical center of the Galaxy. 

The radial distribution of observed and extinction-corrected 1.9 fim flux density is 
shown in Fig. 11. We removed the very bright star IRS 7 and the bins in the distributions 
are independent. The greatest difference between the observed and extinction-corrected 
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distributions occurs at radii outside 10" where the larger extinction corrections associated 
with the CND raise the extinction-corrected flux by a greater amount. The large scale 
near-infrared light distribution was originally measured by Becklin & Neugebauer (1968) 
who fitted an R^'^ power-law at radii between 6" and 10' . The extinction-corrected 
distribution shown in Fig. 11 is clearly somewhat steeper than R^'^ and closer to the Ri-S^O-^ 
power-law shown beside the data. This flux distribution corresponds to a radial density 
distribution p oc R~^-^ in this outer region. This distribution is considerably flatter than 
that derived from number counts of individual stars (oc R~^, Genzel et al. 1996, Allen 
1994, Rieke & Rieke 1994)). These steeper dispersions are strongly influenced by the 'ring' 
of early-type and giant stars at radii of 4-10 " (Genzel et al. 1996). The outer radial flux 
distribution (corrected for extinction) is fit by the power law but a noticeable excess fiux 
is seen in the range 0.8 to 10" radii in Fig. 11. This excess is probably partially due to the 
young stars in the IRS 16 cluster. 

The mean surface brightness distributions of the diffuse starlight are shown in Fig. 12. 
Here, to avoid the strong effects of bright individual stars, the surface brightness (both 
mean and median pixel brightness) were computed from the pixels in the lower 80% of the 
brightness distribution at each radius. (The results are insensitive to the adopted 80% cutoff 
since cutoff between 50 and 85% were tested with similar results; the 80% cutoff was adopted 
finally since it makes use of a large fraction of the pixels. The median also attenuates the 
effects of a few very bright pixels.) Both the observed and extinction-corrected brightness 
distributions increase at smaller radii but exhibit a sharp drop at 0.8" radius, corresponding 
to 0.03 pc or 1.0x10^^ cm. To evaluate whether this drop is statistically significant, we 
independently evaluated the mean and medians for the 4 quadrants (NE, SE, SW, NW) 
from Sgr A* and then estimated the uncertainties from the standard deviation of these 
independent samples - these are the vertical error bars shown in Fig. 12. The central drop 
is significant given these uncertainty estimates. 

The density of discrete stars (as opposed to the total or diffuse light distribution) has 
been derived as a function of projected radius by Genzel et al. (1996). They also found a 
decrease in the projected surface density of star counts. For stars with K < 12 mag, the 
drop occurs at 3" radius; for brighter star with K < 10.5 mag, the decrease occurs inside 
5" radius. Thus our results differ significantly from theirs, showing the drop occuring at ~ 
0.8" . Genzel et al. (1996) attribute the decrease at small radii as being due to coUisions of 
the giant stars with main sequence stars which then deplete the giant stars because of their 
relatively large cross-sections. 

To test whether the drop in the surface brightness distribution inside 1" is due to a 
lack of bright giant stars inside 0.8" , we also show in Fig. 12, the median pixel surface 
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brightness. This median flux should be virtually unaffected by the bright stars which 
dominate a relatively small fraction of the pixels; it still exhibits a modest decrease at small 
radii, although not as large as that seen in the mean surface brightness. 

The drop or levehng off of surface brightness inside 1" could have several explanations: 
it might be the core radius of the nuclear stellar distribution or it might be due to the 

depletion of late type stars with high mass-to-light ratios (Phinney 1989; Bailey & Davies 
1999). Gravitational lensing would produce a drop inside the Einstein radius which is on 
the scale 0.01" (Wardle & Yusef-Zadeh 1992) and is therefore not a likely explanation. 
Sellgren et al. (1990) found a diminishing depth of the 2.3/im CO stellar absorption feature 
in the central few arcsecs, possibly indicating a depletion of giant stars at the center. An 
alternative possibility is ejection of stars from the central 1" by a past or present binary 
black hole (Milosavljevic & Phinney - private communication). 

7. Summciry 

We have used a mosaic image of the 1.87/im HI Pa line to estimate the extinction 
distribution in the Galactic nucleus by comparison with the 6-cm radio continuum and to 
derive the spatial distribution of ionized gas in the central 40" at 0.2" resolution (Nyquist 
sampled in the central 19" ). 

1) In front of the ionized gas, Ay varies from 20 to 50 mag with a median pixel value of 
Ay = 31.1 mag, in excellent agreement with earher estimates derived for the stellar sources. 

2) A broad minimum in the extinction extends ~ 30" NE-SW over the near the 
brightest ionized emission (Northern arm, IRS 16 region, E-W Bar) and increases in the 
extinction are seen along the periphery of the ionized gas towards the dense CND molecular 
features. 

3) The minimum Lyman continuum emission rate from within 40" radius of Sgr A* is 
3.9 xlO^'' sec^^ or Llj/C - 2.7 X 10^ L0. At least half of this emission is within 20" radius. 

4) Mass-loss winds from the hot emission line stars are probably responsible for 
containment of many of the high brightness the ionized filaments. 

5) For both the ionized gas and stellar continuum, the centroids of the large scale 
emission (corrected for extinction derived from Pa) remain within ~ ± 2" of Sgr A*from 
2" out to 40" radius. This provides very strong evidence that Sgr A* (the putative central 
black hole) is indeed at or extremely close (< 0.04 pc or 1.2 x 10^ cm) to the center of the 
Galactic nucleus stellar distribution and presumably the dynamical center of the Galaxy. 
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6) The extinction-corrected enclosed stellar continuum flux is fit quite well by a R^-^ 
power-law at large radii with an excess from 1 to 8" radius. The 1.9/im surface brightness 
increases strongly into 1" radius and exhibits a drop with 0.8" of Sgr A*, possibly 
indicating the influence of the central black hole on the nearby stellar distribution. 

This paper is based on observations with the NASA/ESA Hubble Space Telescope 
obtained at the Space Telescope Science Institute, which is operated by Association of 
Universities for Research in Astronomy, Incorporated, under NASA contract NAS5-26555. 
The NICMOS project was supported by NASA grant NAG 5-3042 to the NICMOS 
instrument definition team. We thank Doug Roberts for use of his H92q; image. It is a 
pleasure to thank Zara Scoville for assistance in editing. 

8. Appendix- Stelleir Continuum Subtraction 

Subtraction of the stellar continuum emission from the 1.87 fim image was one of the 
most challenging aspects of this study. Compounding the problem of continuum subtraction 
from the hne-|-continuum image (F187N), the stellar sources have different intrinsic colors 
and variable reddening. In addition, some stars have Pa emission or absorption. 

For subtraction of the stellar continuum from the F187N image, we experimented with 
three techniques: fixed-scale factor subtraction of F190N from F187N; PSF fitting and 
subtraction of the stars from the F187N image; spatially variable scale factor subtraction of 
F190N from F187N. We discuss these approaches separately in the following subsections. 

8.1. F187N - CxF190 

The most straightforward approach is to multiply the F190N image by a fixed scale 
factor and subtract this from the F187N image, with the scale factor being that which 
yields the smallest positive and negative residuals at the locations of identified stars. Prom 
trials with the scale factor varying between 0.85 and 1.0 in steps of 0.01, we determined 
that the best overall scale factor was 0.91 and the resultant residual image is shown in 
Pig. 13. This scale factor is confirmed in the peak of the histogram of F187N/F190N ratios 
shown in Fig. 14. The bright positive and negative spots in the residual image are due to 
too little or too much stellar continuum being subtracted for the individual stars. Positive 
residuals predominate in the central region, negative residuals near the periphery. This 
systematic variation clearly indicates that the reddening is lowest in front of the central 
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area and increases over the periphery. Thus a straightforward subtraction of the F190N 
image (niTiltiphed by a fixed-scale factor) from the F187N image leaves strong negative 
and positive residuals at the positions of the stars. This implies that suitable continuum 
subtraction requires different scale factors at the location of each individual star. 

8.2. Stellar PSF Subtraction 

The second approach removed the stars directly via PSF subtraction. For the PSF 
fitting we employed the 'starfinder' IDL software package (Diolaiti et al. 2000). The 
software is designed for imaging with oversampled PSFs; thus, PSF fitting was only 
attempted for the NIC2 mosaic, which only marginally meets this criterion. A composite 
PSF was derived from the NIC2 mosaic itself and was dominated by the brightest source 
in the field, IRS7. PSF subtraction was attempted both by subtracting an image of the 
detected stars (detected fluxes convolved with the PSF) from the F187N directly as well 
as subtracting stars detected at both F187N and F190N from the fixed scale factor Pace 
image. We found that PSF fitting with 'starfinder' worked satisfactorily for deriving stellar 
fluxes in the F187N and F190N images which was used to construct a list Pa emission 
stars. However, we found that the PSF subtraction of the reconstructed point sources did 
not give as satisfactory results in removing stellar residuals, possibly due to a non-ideal 
PSF derivation coupled with the only marginal sampling. This undersampling means that 
pixelation errors which are quite unpredictable become important. In addition, the strongly 
varying diffuse Pa background over both small and large scales makes the PSF fitting less 
accurate than if the field were entirely composed of point sources. We therefore decided 
after much effort that PSF fitting and subtraction would not yield acceptable results in 
removing the stellar continuum. 

8.3. F187N - C(x,y)xF190 

Lastly, we experimented extensively with a procedure which determined individual 
scale factors for approximately 1000 stars identified in the F190N and F187N images, used 
these scale local factors in the vicinity of each of the stars and then adopted the average 
stellar scale factor for all remaining pixels. This produced a final residual image, shown in 
Fig. 2 considerably cleaner than that obtained using a fixed scale factor (Fig. 13), although 
the procedure involved a number of complex steps. 

Specifically, : 
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1) On both the F187N and F190N images, we first removed the diffuse low-level 
background (Pa and faint stars) determined as the median intensity in the 16 x 16 pixel 
area centered on each original 0.0753" pixel. (The value of 16 pixels was chosen to be 
enough that at least 50% of the pixels would be outside the PSF of bright stars, yet small 
enough that variations in the background could be picked up.) 

2) The F190N pixels associated with each stellar peak were then collected using the 
algorithm described by Scoville et al. 2001. In order to be considered as a star, the peak 
pixel had to exceed 410 //Jy in the F190N image; and the adjacent pixels were included in 
the star down to 175 /xJy. Thus the pixel boundaries of over 1500 stars were determined. 

3) Local scale factors for each star were then derived as the ratio of the total flux 
measured in the F187N and F190N images within identical pixel boundaries for each star 
defined from F190N. The derived scale factors for 1031 stars are shown in Fig. 14. 

4) An image of the continuum subtraction scale factors (C(x,y)) was then made 
by setting the local scale factor within each star's pixel boundaries to be that derived 
individually in step 3, and setting the scale factor for all pixels not included in any star 
to the average of all the stellar scale factors.^ The use of the 'average' scale factor in the 
'non-stellar' pixels does not introduce serious errors - the precise scale factor there makes 
little difference since the amount of F190N signal which is to be subtracted is lower than 
the Palpha in F187N in these areas. 

5) A difference image was then constructed as F187N - C(x,y) xF190N. 

6) The difference image was generally better in terms of lower residuals (compared to 
those in Fig. 13), but still had small residuals near some of the brighter stars (due probably 
to pixelation errors). Therefore, above the 410 /^Jy in the stellar positions, we interpolated 
from the immediately surrounding pixels.^ 

Although the above process was complex and the last step (# 6) is obviously ad 
hoc, the resulting image retains the maximum areal coverage, replacing only the brightest 
stellar pixels with interpolated values, and results in an image relatively free of residual 
stellar effects. The focus of the present paper is the interstellar ionized gas and extinction 
distributions. The ionized gas in the directions of the stars is clearly more accurately 



^Any stars with apparent F187N/F190N flux ratios greater than 1.1 -a 20% effect since the peak scale 
factor is 0.91-were deemed to probably be strong Pa emission line stars (see below); for these stars, the local 
scale factor was also set to the average scale factor of the remaining stars. 

^In the vicinity of the very bright star IRS 7 which had a peak continuum surface brightness ^ 30,000 
/uJy per pixel, the area with strong PSF efl'ects (30 x 25 pixels) was entirely removed and interpolated. 



studied with high resolution spectroscopy. 
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Fig. 1. — The mosaiced F187N and F190N images are shown with pixels sampled to 0.0753" 
and rotated with north up. The region is dominated by stellar sources, and the Pa emission 
is seen as relatively faint diffuse emission in the F187N image. The color scale units are 
/xJy per 0.753" pixel and the offsets are relative to Sgr A*(a2ooo = 17''45™40'*.058, 62000 
= -29° 00'27.90"). In the f.87;um image, contours of HCN (1-0) emission are shown and 
the boxes indicate the positions of the Pa emission stars listed in Table 2. The HCN 
countours are dashed when negative (i.e. absorption) and the levels are at : -100,-50,-25,- 
10,-5,2,4,6,8,10,15,20,25,30,35,40,45,50,60 Jy km/ sec /beam. 

Fig. 2. — Pa emission for the entire region imaged. The color scale units are 10~^^ ergs cm~^ 
sec~^ per 0.0753" pixel and spatial offsets are relative to Sgr A*. 

Fig. 3. — The extinction distribution derived from the observed ratio of Pa to 6-cm radio 
continuum emission. The Pa emission was convolved to the same spatial resolution as 
the radio and we assume the radio continuum is entirely free-free except on Sgr A*. The 
contours (-100,-50,-25,-10,-5,2,4,6,8,10,15,20,25,30,35,40,45,50, and 60 Jy km s'^) are for the 
HCN emission at 3" resolution (Christopher et al. 2003). A broad minimum in the extinction 
is seen in the IRS 16 cluster and the extinctions increase at the periphery of the ionized gas 
in the molecular features. The dashed contours in the HCN indicate locations where the 
HCN appears in absorption of the radio free-free (from the ionized gas) and nonthermal (Sgr 
A*) continuum. The apparently high extinction on Sgr A* is due to the strong non-thermal 
contribution which yields spuriously large extinction estimates. Likewise, the interstellar 
extinction towards the bright Pa emission line stars cannot be accurately measured; the 
excess photospheric Pa emission yields spuriously low extinction values towards these stars. 

Fig. 4. — Histogram of visual extinctions (Ay) for pixels with detected radio continuum and 
Pa emission. The median extinction is Ay = 31.1 mag. 

Fig. 5. — The extinction distribution derived from the observed ratio of Pa to H92a radio 
recombination hne emission (Roberts & Goss 1993). The Pa emission was convolved to the 
same spatial resolution as the radio (2" ). The contours are for the HCN emission at 3" 
resolution (Christopher et al. 2003). This extinction distribution agrees closely with that 
derived independently from the 6-cm radio continuum (Fig. 3). A broad minimum in the 
extinction is seen in the E-W bar and Northern arm and the extinctions increase at the 
periphery of the ionized gas in the molecular features. We note that regions with insufficient 
signal/noise in the H92a image are masked out (such as the region to the SW of Sgr A*). 
The dashed contours in the HCN indicate locations where the HCN appears in absorption 
of the radio free-free (from the ionized gas) and nonthermal (Sgr A*) continuum. 
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Fig. 6. — Histogram of the ratios of extinctions determined from Pa/E.92a versus those 
determined from Pa/6-cm radio continuum. The histogram is of course restricted to those 
pixels shown in both Figs. 5 and 3 with determinant extinctions. 

Fig. 7. — The Pa emission corrected for extinction using the extinction derived from the 
ratio of Pa to radio continuum emission (Fig. 3). The color scale units are 10~^^ ergs cm~^ 
sec~^ per 0.0753" pixel and spatial offsets are relative to Sgr A*. The contours are for the 
HCN emission at 3" resolution (Christopher et al. 2003). 

Fig. 8. — Histogram of extinction corrected Pa fluxes (per 0.0753" pixel) for pixels with 
extinction estimates derived from the ratio of Pa emission to radio continuum. 

Fig. 9. — The Pa intensity distribution along an E-W strip (1 pixel wide) crossing the 
Northern arm is shown. The peak occurs at 7.2"E, 11"N of Sgr A*. The two traces are for 
the observed (solid line) and extinction-corrected (dashed line) Pa fluxes. 

Fig. 10. — Position centroids of extinction-corrected Pa (small boxes) and extinction- 
corrected 1.9 /im (large boxes) emissions are plotted on the observed 1.9 /xm emission within 
± 2.5" of Sgr A*. The centroids were measured successively for pixels within 2, 5, 10, 15, 25, 
30, 35, 40" radii of Sgr A* and each box is labelled with the corresponding radius. The very 
bright source IRS 7 was excluded and the extinction within ± 0.5" of Sgr A* was set equal 
to the average in the annulus between 0.5 and 0.75" radius (on account of the non-thermal 
flux of Sgr A* which precluded extinction estimates directly on Sgr A*). 

Fig. 11. — The radial distribution of observed and extinction-corrected 1.9 /xm flux densities 
are shown as a function of aperture radius. The very bright source IRS 7 was excluded and 
the extinction within ± 0.5" of Sgr A* was set equal to the average between 0.5 and 0.75" 
radius (on account of the non-thermal flux of Sgr A* which precluded extinction estimates 
directly on Sgr A*). Also shown are an R^-^ power-law which Becklin & Neugebauer (1968) 
flt to the observed 2.2 //m flux distribution on larger angular scales (6" to 10' ). The R^-^ 
power-law flts somewhat better to the extinction-corrected 1.9 //m flux distribution. 

Fig. 12. — The surface brightness distribution for the observed and extinction-corrected 1.9 
//m radiation is shown. Discrete, bright sources (such as IRS 7) were excluded by removing 
the brightest 20% of the pixels in each radial bin. Because the noise is very non-Gaussian 
(dominated by the stellar brightness distribution), the vertical error bars were estimated 
from the standard deviation of values calculated in the 4 quadrants relative to Sgr A*. 

Fig. 13. — The image of Pa emission obtained by subtraction of 0.91 xF190N from the F187N 
image. The bright positive and negative spots in the residual image are due to too little 
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or too much stellar continuum being subtracted for the individual stars. Note that positive 
residuals are prevalent in the central regions and negative residuals near the periphery - a 
result of the fact that the reddening is lowest in front of the central area and increases in 
the periphery. In addition, the Pa emission line stars (see Table 2) are concentrated in the 
central region. 

Fig. 14.— The distribution of measured stellar flux ratios F187N/F190N is shown for 1031 
stars determined with the spatially variable scale factor method described in the Appendix. 
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Table 1. Galactic Center Pa Observations 



Camera 


Targef^ 


Date 


Filter 


A (/im) 


Exp. Time (s) 


Orient, angle 


NICM0S3 


GC-A 


Apr 6, 1998 


F187N 


1.874 


40 X 4^ 


24.54 


NICM0S3 


GC-A 


Apr 6, 1998 


F190N 


1.900 


40 X 4 


24.54 


NICM0S3 


GC-B 


Apr 6, 1998 


F187N 


1.874 


40 X 4 


24.54 


NICM0S3 


GC-B 


Apr 6, 1998 


F190N 


1.900 


40 X 4 


24.54 


NICM0S3 


GC-C 


Apr 6, 1998 


F187N 


1.874 


40 X 4 


24.86 


NICM0S3 


GC-C 


Apr 6, 1998 


F190N 


1.900 


40 X 4 


24.86 


NICM0S3 


GC-D 


Apr 6, 1998 


F187N 


1.874 


40 X 4 


24.86 


NICM0S3 


GC-D 


Apr 6, 1998 


F190N 


1.900 


40 X 4 


24.86 


NICM0S2 


GC (center) 


Aug 29, 1998 


F187N 


1.874 


304 X 4 


223.81 


NICM0S2 


GC (center) 


Aug 29, 1998 


F190N 


1.900 


304 X 4 


223.81 



^All data from Program ID 7250 

^Four dithers (spaced 1.1" apart in a square pattern) executed for each target 

^J2000 positions of field centers: GC-A(17:45:41. 9,-29:00:04), GC-B(17:45:41.9,-29:00:52),GC- 
C(17:45:38. 1,-29:00:52), GC-D(17:45:38.1,-29:00:04),GC-Center(17:45:40.0,-29:00:28) 

■^Position angle E of N 



TABLE 2 

Pq Emission Line Stars in the Central Parsec'* 



Source 


A« (")'' 


AS (")'' 


Ar (")'' 


F187N/F190N ^ 


d 


Fp„ ■= 


Identification * 


1 


1.13 





50 


1 


24 


1 


36 


10.9 


215 


IRS 16C 


2 


0.00 


1 


24 


1 


24 


1 


20 


11.0 


107 


IRS 16NW 


3 


1.00 


-0 


94 


1 


38 


1 


41 


11.0 


216 


IRS 16SW 


4 


-0.97 


2 


09 


2 


31 


1 


10 


13.3 


7 




5 


2.86 


1 


10 


3 


06 


1 


23 


9.9 


329 


IRS 16NE 


6 


2.91 


-1 


16 


3 


13 


2 


59 


13.1 


126 


IRS 16SE2 


7 


0.61 


-3 


09 


3 


15 


1 


21 


11.0 


106 


IRS 33SE 


8 


-3.20 


-1 


43 


3 


50 


2 


15 


12.7 


122 


IRS 13E4 


9 


-3.15 


-1 


72 


3 


59 


2 


07 


11.9 


236 


IRS 13E2 


10 


-4.05 


1 


62 


4 


36 


1 


36 


13.7 


15 


IRS 34 


11 


2.92 


3 


51 


4 


56 


1 


62 


14.4 


14 


HeIN3 


12 


-4.37 


-1 


64 


4 


67 


1 


53 


13.6 


25 




13 


-3.71 


2 


87 


4 


70 


1 


36 


14.3 


9 




14 


0.91 


5 


81 


5 


88 


1 


18 


12.7 


19 




15 


2.81 


-5 


59 


6 


26 


2 


07 


13.0 


89 


IRS 9W 


16 


-3.95 


4 


95 


6 


33 


1 


63 


13.1 


47 


IRS 7W 


17 


-2.46 


6 


01 


6 


49 


1 


80 


14.3 


19 




18 


4.41 


5 


02 


6 


69 


1 


84 


13.8 


32 


IRS 7E2 


19 


-7.63 


-3 


57 


8 


42 


1 


88 


12.5 


113 


AFNW 


20 


4.22 


-7 


97 


9 


02 


1 


77 


14.2 


21 




21 


9.47 





37 


9 


48 


1 


92 


13.6 


43 


ID180 


22 


-6.56 


-6 


92 


9 


54 


3 


09 


11.0 


479 


AF 


23 


-9.64 


-2 


57 


9 


98 


2 


38 


13.3 


90 




24 


-8.64 


-5 


36 


10 


17 


1 


75 


13.8 


32 




25 


-1.53 


10 


05 


10 


17 


1 


78 


12.7 


81 


IRS 15SW 


26 


1.24 


11 


81 


11 


88 


1 


38 


12.0 


80 


IRS 15NNE 



^Stars for which the flux in F187N is > 10% more than in F190N (where the F190N image was multiplied by a factor 
of 0.91 as discussed in the text) 

"^Offsets and radius are relative to Sgr A* (^2000 ^ 17'^45"^40^ .058, ^9000 ^ -29° 00'27.90"). 
'^FISTN/FIQON flux ratio 

'^FIOON flux density (magnitude), not corrected for extinction; zero magnitude — 808 Jy. All fluxes extracted from 
NIC2 images using "starflnder" with the exceptions of IRS 13E4 and IRS 13E2 (which required PSF subtraction) and 
IRS 15NNE (which appears only in the NIC3 image, and for which aperture photometry was performed). 

'^Units : 10~^^ ergs cm~^ sec~^. The fluxes have not been corrected for extinction. 

^Identifications with stars listed by Genzel ei al. (2000), Ott, Eckart &z Genzel (1999), Paumard et al. (2001), and 
Maillard et al. (2003). Our positions were checked against those given by Mentcn et al. (1997) and we find agreement 
to within 0.15 in all cases. 
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